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We study the cosmology of Galileon modified gravity models in the linear perturbation regime. We de- 
rive the fully covariant and gauge invariant perturbed field equations using two different methods, which give 
consistent results, and solve them using a modified version of the CAMB code. We find that, in addition to mod- 
ifying the background expansion history and therefore shifting the positions of the acoustic peaks in the cosmic 
microwave background (CMB) power spectrum, the Galileon field can cluster strongly from early times, and 
causes the Weyl gravitational potential to grow, rather than decay, at late times. This leaves clear signatures 
in the low-/ CMB power spectrum through the modified integrated Sachs-Wolfe effect, strongly enhances the 
linear growth of matter density perturbations and makes distinctive predictions for other cosmological signals 
such as weak lensing and the power spectrum of density fluctuations. The quasi-static approximation is shown 
to work quite well from small to the near-horizon scales. We demonstrate that Galileon models display a rich 
phenomenology due to the large parameter space and the sensitive dependence of the model predictions on the 
Galileon parameters. Our results show that some Galileon models are already ruled out by present data and 
that future higher significance galaxy clustering, ISW and lensing measurements will place strong constraints 
on Galileon gravity. 



I. INTRODUCTION 

The accumulated evidence for the present-day accelerated 
expansion of the Universe, driven by what is generically re- 
ferred to as 'dark energy', is now overwhelming JH-d]. The 
simplest explanation for the nature of dark energy is a simple 
cosmological constant but, despite the good agreement with 
the observational data so far, such an explanation is plagued 
with serious fine tuning and coincidence problems. This has 
motivated the proposal of alternative models to explain the 
observations, the majority of which fall into two classes. The 
first one assumes the existence of a dynamical dark energy 
field (often of scalar type) which dominates the energy density 
today and has a negative pressure to accelerate the Universe 
0, Si . The other considers that the standard law of gravity, 
general relativity, fails on cosmological scales and must be 
completed by modifications capable of accelerating the Uni- 
verse J(|. Models in the second class have attracted a lot of 
research interest recently, and significant progress has been 
made in both the theoretical modelling |0, H and numerical 
simulations HHU]. 

One notable example of a modified gravity model which 
has been the subject of many recent papers is the Galileon 
model lfl2l [kill . Here, the deviation from general relativity 
is mediated by a scalar field <p, dubbed the Galileon, whose 
Lagrangian is invariant under the Galilean shift symmetry 
d^cp — > d^ip + (hence the name), where is a constant 
vector. Such a field appears, for instance, as a brane-bending 
mode in the decoupling limit of the DGP braneworld model 
Il4l - [l6ll which was proposed well before the Galileon model. 
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However, in spite of being theoretically appealing, the self- 
accelerating branch of the DGP model, which is of interest to 
the cosmological community, is plagued by the ghost prob- 
lem 1IT7U19I1 (i.e. there is not a well defined minimum en- 
ergy). It was shown that by generalising the DGP Lagrangian, 
in the decoupling limit, to more general forms that satisfy 
the Galileon shift symmetry it is possible to avoid the pres- 
ence of ghosts Il2tl . In a flat Minkowksi spacetime there are 
only five possibilities for the Galileon Lagrangian. These in- 
clude highly nonlinear derivative couplings of the field but that 
lead only to second-order field equations due to the protection 
offered by the shift symmetry. In 11311 . a generalisation of 
these Lagrangians to curved spacetime was presented by in- 
troducing explicit couplings between the Galileon field and 
the spacetime curvature. Such couplings change the way in 
which spacetime curvature responds to matter distributions, 
which is why the Galileon model is a subclass of modified 
gravity theories. 

In any viable modified gravity theory, it is crucial that de- 
viations from standard gravity get suppressed (or screened) 
in high matter-density regions where general relativity has 
been tested to high accuracy rf20l, I2II1 . In the case of Galileon 
gravity, such a screening is realised via the Vainshtein mech- 
anism Il22ll . which relies on the presence of the nonlinear 
derivative self-couplings of the Galileon field. Here, far away 
from gravitational sources, the nonlinear terms are subdom- 
inant and the Galileon field satisfies a linear Poisson equa- 
tion (as the Newtonian potential), so that the extra (fifth) force 
mediated by it can be sizeable and proportional to standard 
gravity, effectively renormalising Newton's constant. Near 
the sources, on the other hand, the nonlinear terms become 
important, which strongly suppress the spatial variations of 
the Galileon field compared to that of the Newtonian poten- 
tial and ensure that the extra force, which is the gradient of 
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the Galileon field, is not felt on scales smaller than a given 
'Vainshtein radius'. In certain respects, this is very similar to 
the chameleon screening 11231 12411 . which operates for instance 
in f(R) gravity models 125142711 . However, in the chameleon 
case the self-interaction of the scalar field depends on the field 
value (through a nonlinear interaction potential) rather than its 
derivatives, and the non-derivative coupling of the scalar field 
to matter makes its behaviour highly sensitive to the environ- 
mental matter density - in high density regions the field value, 
rather than merely its gradient, becomes extremely small so 
that the extra force is suppressed. 

It is therefore evident that one has to go beyond the local 
environment to look for possible deviations from general rel- 
ativity and distinct signatures of the different modified grav- 
ity models. In particular, a promising way is to look at the 
cosmic expansion and the formation of structure in the Uni- 
verse: different screening mechanisms in different modified 
gravity models can lead to very different predictions as to 
when, where and how the various cosmological observables 
are affected. 

The effects of Galileon gravity models on the background 
cosmological expansion have already been studied in the lit- 
erature in great detail I28143311 . It has been shown that in these 
models there is a stable de Sitter point that can be reached 
after the radiation and matter dominated eras, thus yielding a 
viable cosmological expansion history. Conditions to avoid 
the ghosts and other theoretical instabilities have also been 
derived by considering the linear perturbations IT30l[32Tl . 

To improve our understanding of the cosmological effects 
of Galileon gravity models and make direct comparisons with 
observational data, a proper investigation of the evolution of 
density fluctuations and formation of large-scale structure is 
necessary. Here, as an initial step, we consider the regime in 
which the density fluctuations are small such that their evo- 
lution is well described by linear perturbation theory. This 
regime is relevant for several important cosmological observ- 
ables, such as the power spectrum of the cosmic microwave 
background (CMB) temperature fluctuations and its polarisa- 
tions, the growth of matter density perturbations, the weak 
gravitational lensing of distant galaxies and the CMB map, 
and the integrated Sachs-Wolfe (ISW) effect and its cross cor- 
relation with the galaxy distribution. The rich information 
contained in this regime therefore warrants a detailed study of 
the Galileon effects, which is precisely the topic of this paper. 
The nonlinear regime of structure formation can in principle 
contain further interesting information, but its study is beyond 
the scope of the current paper. 

The layout of this paper is as follows. We start by briefly 
presenting the Galileon model and the Galileon and metric 
field equations of motion in Section[II] The perturbation equa- 
tions are derived and presented in a covariant and gauge in- 
variant (CGI) way in Section [TTT] using the method of 3 + 1 
decomposition. In Appendix IE1 we present an alternative and 
considerably simpler derivation of the perturbation equations 
which is particularly suitable for the Galileon model as it takes 
advantage of the fact that the Lagrangian density is fixed by 
the Galilean shift invariance and that there are no derivatives 
higher than second order. We present and discuss the results 



for the CMB, lensing and linear matter power spectra in Sec- 
tion [TV] which we obtain using a version of the CAMB code 
J34ll which we have modified. In Section [TV] we also discuss 
the time evolution of the gravitational potential, Galileon field 
perturbation and Galileon density contrast and the validity of 
the quasi-static limit. We conclude in SectionlVl 

Throughout this paper we will use the unit c = 1 and metric 
convention (+,—,—,—). Greek indices run over 0,1,2,3 and 
we will use 8ttG = k = Mpj interchangeably, where G is 
Newton's constant and Mpi is the reduced Planck mass. 



II. THE MODEL 



The covariant uncoupled Galileon action can be written as 



S 



R 



16nG 
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(l) 



where g is the determinant of the metric, R is the Ricci 
scalar and ci_5 are dimensionless constants. The five covari- 
ant terms in the Lagrangian densities, which are fixed by the 
Galilean invariance, d^ip —> d^ip + 6,,, are given by 

£1 = M 3 v?, 

£2 = v^VV, 

U = -^V^V> [2(L» 2 - 2(V M V^)(V J V» 

£5 = -^V^V> [(L» 3 - 3(D V )(V M V^)(VV» 

+2(V M V»(V t ,V'V)(V p V'V) 
-6(V^)(VV»(V^)G„,] , (2) 

where ip is the Galileon field and M 3 = Mp\Hq with Hq 
being the present-day Hubble expansion rate. 

Although the form of the Galileon Lagrangians, £;, is fixed 
by imposing the Galilean invariance, there is still freedom 
to introduce a coupling between the Galileon and the matter 
fields. For instance, 113211 considers a derivative coupling in the 
Einstein frame of the form C coup u ng ~ T^V uipV v (p, which 
translates into C coup ii ng ~ C 1 " V M </V ' v <p in the Jordan frame 
(see e.g. I3~5l437ll ). where T^ u and C" are the matter energy- 
momentum tensor and the Einstein tensor, respectively. 

In the rest of the paper we will choose to work in the Jordan 
frame adding to Eq. (Q} the Lagrangian density 



£g = -cg^G^V^V^, 



(3) 



where cq is a dimensionless constant which determines the 
strength of the coupling. We will be interested in the cases 
where the acceleration is due only to the field kinetic terms 
and therefore we will set the potential term c\ to zero. 
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The modified Einstein equations and the Galileon equation 
of motion are obtained by varying the action, S, with respect 
to g^y and <p, respectively. Our derivation agrees with those 
present in the literature lfl3l l32ll although we explicitly write 
the Riemann tensor in terms of the Ricci and Weyl tensors, 
which leads to the cancellation of some terms and hence to a 
slight simplification of the final expressions. We show these 
equations in Appendix lAl 



III. THE PERTURBATION EQUATIONS 
A. The Perturbed Equations in General Relativity 

In this section we derive the covariant and gauge invariant 
perturbation equations in Galileon gravity. This will be done 
in detail below but before that let us outline the main ingredi- 
ents of 3 + 1 decomposition and their application to general 
relativity for ease of later reference. 

The main idea of 3 + 1 decomposition is to make spacetime 
splits of physical quantities with respect to the 4-velocity w M 
of an observer. The projection tensor is defined by — 
9nv — u^u v and can be used to obtain covariant tensors which 
live in 3-dimensional hyperspaces perpendicular to vf 1 . For 
example, the covariant spatial derivative V of a tensor field 



T 



/3-T 



is defined as 



v°t£7 = h a X ■ ■ ■ KK ■ ■ ■ 1 (4) 

The energy-momentum tensor and covariant derivative of 
the 4-velocity are decomposed, respectively, as 



Tpv = ir^ + 2q (jJi u v ) + pu^u v - ph^, (5) 
V^Uv = o + t^^u + T^Oh^, + u^A Ul (6) 

where 7r„j, is the projected symmetric and trace-free (PSTF) 
anisotropic stress, is the heat flux vector, p is the isotropic 
pressure, the PSTF shear tensor, vo^ = V^u^ the vor- 
ticity, 9 = V Q M a = 3d/ a = 3H (a is the mean expansion 
scale factor) the expansion scalar and = ii^; the overdot 
denotes a time derivative expressed as </> = u a W a <f>, brack- 
ets mean antisymmetrization and parentheses symmetrization. 
The normalization is such that u a u a = 1, The quantities ir^, 
q^, p and p are referred to as dynamical quantities and cr pl/ , 
tu^ u , 9 and as kinematical quantities. Note that the dy- 
namical quantities can be obtained from Eq. (0 using the re- 
lations 



P=--h^ u , 



h"u p T u 



7r M „ = h^hlTpr +ph llv . 



(7) 



Decomposing the Riemann tensor and making use of Ein- 
stein equations, we obtain, after linearization, five constraint 
equations J38ll : 



= V 



a ( 



vrvo. 



2V, 



B 7 

e \ u . 



(8) 

(9) 
(10) 

(11) 



V V B^ = i« 



V a q,3 + (p+p)^a/3 



^ p W; (12) 



and five propagation equations: 

= 



±e 2 -V-A+^(p + 3 P ), 



2 - K 



= - 



(p + p)(X^ v +V {tl q v) 



(13) 
(14) 
(15) 

(16) 
(17) 



Here, e MI/Q( g is the covariant permutation tensor, £^ v and 
B^y are, respectively, the electric and magnetic parts of the 
Weyl tensor W M „ Q/ 9, defined by E^ v = u a «% (m ^ and 
B^v = —^"u^e^jfWjSvfi- The angle brackets mean tak- 
ing the trace-free part of a quantity and V • v = V a v a , where 
v is an arbitrary vector. 

Besides the above equations, it is useful to express the pro- 
jected Ricci scalar R into the hypersurfaces orthogonal to w M 
as 



R = 2np- ^9 2 . 



(18) 



The spatial derivative of the projected Ricci scalar, rj^ 
aV uR/2, is given as 



2a - 

j]^ = naV^p - — QV^O, 
and its propagation equation given by 



(19) 



29 



2a9 



-V^V-A-aKV^V-g. (20) 



Finally, there are the conservation equations for the energy- 
momentum tensor: 

p+(p + p)9 + V -q = 0, (21) 
% + ^% + + P)\ - V M p + V 1 ^ = 0. (22) 
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In this paper we will always consider the case of a spatially- 
flat Universe and, as a result, on large scales the spatial curva- 
ture vanishes R = 0. Thus, from Eq. ( fTSI l. we obtain the first 
Friedmann equation 



Kp. 



(23) 



Note that at the background level only the zeroth-order terms 
contribute to the equations. The second Friedmann equation 
and the energy-conservation equation are obtained by taking 
the zeroth-order parts of Eqs. ( [T3ll2T1 i. as 



the Ricci tensor in terms of the kinematical quantities. 
For this let us expand the symmetric rank-2 tensor R^ in the 
following general way 



R 



(26) 



in which T M is a four-vector and a PSTF rank-2 tensor, 
both of which live in the 3-dimensional hyperspace perpendic- 
ular to the observer's four-velocity (u M T M = « M £ M „ = 0). A 
and H are scalar quantities. Then, using the modified Einstein 
field equations 



(27) 



\o 2 + ^(p + 3p) =0, 
p+(p + p)6 = 0. 



one gets, 



(24) 
(25) 



In what follows, we will only consider scalar modes of per- 
turbations, for which the vorticity, w^ u , and the magnetic part 
of the Weyl tensor, £> M „, are at most of second order l38ll and 
will be neglected from our first-order study. 



B. The Perturbation Quantities in Galileon Gravity 

In the effective energy-momentum tensor approach, the 
field equations Eqs. d8l-[25Tl above preserve their forms, but 
the dynamical quantities p, p, g„ and ttuv should be replaced 

/ 



by the effective total ones p tot = p' 



1u 



ql + q G and tt^ 4 



\1,V 



p G , p tot 



P J 



P 



, in which the super- 



scripts G and ' identify the contributions from the Galileon 
field and the rest of the matter fluid (including cold dark mat- 
ter, baryons, photons and neutrinos), respectively. From here 
on we shall drop the superscript tot for ease of notation. 

Before using Eq. (0 to calculate p G , p G , q G and ir G 
from the components of the Galileon energy-momentum ten- 
sor shown in Appendix[A], we need an explicit expression for 



A = 



T, 



-n(p + 3p) 



-7; K (P-P) 
1 

k<7m 



V • A 



R-V A 



2V„0 



2 



///,' ) 



V <A^> 



(28) 



(29) 



(30) 



(31) 



where we have used Eqs. (f9l IT3l IT4l IT8T >. Notice that the first 
lines are expressed in terms of total dynamical quantities and 
the second lines in terms of kinematical quantities. 

With the above useful relations and after some tedious but 
straightforward calculations, the Galileon contribution to the 
energy-momentum tensor can be identified as 



<-2 



C3_ 

M 3 



2p 3 9 + 2ip 2 Up 



M 6 



^e 2 + Ap 3 eup + ^r 



C5 

M 9 

ri 



, Mpi 



^ 2 d 2 + ^en^ + ^ 2 R 



<-2 



2* 



C4 

M 6 

_C5_ 

M 9 
Mpi 
' M 3 



-A(pLp 3 e - ip 4 e - i^> 4 (9 2 - A(ptp 2 ntp - -<p 3 6n<p + <^ 4 v ■ a + ^v? 4 -R 



5 •• -Art 2 ,-J>i 



Cg 



-(pip 

4 

~3 



1 - 5q3 2 ip i e 2 nip - ->pip 3 e\jip - -p^R - -^enp + -ifev ■ a 



9 



<P 



9' 



C2 



ipV^ip 



_C3_ 

M 3 



hp 2 e - hp 2 e 2 + ^p 2 v . a - ^D(p - ^(peup + ^(p 2 R 



2<p 2 6V f t<p - 2<£ 2 V„0 



(32) 



(33) 
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Ca_ 

M 6 

C5 

M 9 



-4ifP6$p<p + 2p 3 9 2 V^<p - ffipO + -^V"^ + -^ 4 V a ^ a 



, Mpi 



C5 

'M 9 



- (ipH + ^e 2 + scptfe) - {^e + 3<^ 4 ) o> - Ucptfe + + \^o 2 ^} v <M v v) ^ 



M- 



pi 



M 3 



c G 



(34) 



(35) 



r 



in which □ = V M V, 



Following the same procure, the perturbed Galileon field 
equation of motion (see AppendixlAl is given by 



= c 2 



M 6 

_C5_ 

'm 9 



_C3_ 

M 3 



ktptpQ + -tpODtp + itpOip + 2tp 2 9 2 + 2if 2 9 - 2(f 2 V ■ A 
3 



26 - 

%(p<p 2 9 2 + itp 3 66 + 2ip 3 6 3 + %<f(peU(p + —p 2 9 2 U<p - 4<y9 3 (9V ■ A + 4(p 2 9Dip + 3</30 2 i? + -<p a 9R 

9 o 



20 
"9 



■ <ptp- 



3/j3 



-0 4 ^ 2 + -<p 3 6> 3 D<^ + -tp 4 9R 
i 9 2 



-<p 4 9 2 R - -o? 4 2 V • A + ^Lp 2 e 2 Uip + -ip 3 99D(p + 2<pip 3 9R 
6 3 3 



M- 



pi 



M 3 



CG 



2 9 4 • - 2 , - 4 
-^ 2 + -^+-0 2 D^+- 



2 , 4 1 
-<^>6> 3 - -y>0V • A + ipR + -tp9R 

3 3 3 



(36) 



As a consistency test, we checked that Eqs. 
the conservation Eqs. (f2Tl [22t . 



satisfy 



^ = E 73^' v = -E^' 



(37) 



C. Perturbed Equations in fc-space 



in which Q k is the eigenfunction of the comoving spatial 
Laplacian a 2 D satisfying 



For the purpose of the numerical studies presented in this 
paper, we need to write the perturbed quantities derived in the 
last subsection in terms of fc-space variables. This is achieved 
with the aid of the following harmonic definitions: 



v*<p = E 

fe 


. L.2 

fe 


A n = E a A ^' 

k 


v,p-E^ 




z — ' a 



(38) 



and Q\ and Q*„ are given by Q\ = f V M Q fc and by Q*„ = 
f V^Q^, respectively. 

In terms of these harmonic expansion variables, Eqs. (l9l[Tl~l 
[l4l[T6l[T9l|2Qb can be rewritten as 



;k 2 (a — Z) = nqa 2 



-no 2 [fc(II + x) + 'iUq\. 



(39) 
(40) 
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k(<7/ + Ha) = k 2 ((f) + A) - -na 2 U : (41) 

k 2 (4>' + H<f>) = ^na 2 [k{p + p)<r + kq-ILf- HU] ,(42) 
k 2 7] = K X a 2 - 2kUZ, (43) 



krjf = -nqa 2 - 2kHA, (44) 

respectively, where % = a 1 /a and a prime denotes a deriva- 
tive with respect to conformal time r (adr = dt, with t the 
physical time). From Eqs. d32l[34l[33T l one obtains the fc-space 
variables x G , 1 G and I1 G 



c 2 l (v,V + ^) + ^ ([18^7' + 18^' 3 HA] + k [2^Z] + k 2 [2^]) 



[90<f' 3 n 2 j' + 90<p /4 H 2 A] + k [l5ip /4 HZ] + fc 2 



C4 1 

([105^ /4 W 3 7 ' + 105^ ,5 W 3 ^] + fc [2V 5 H 2 Z] + fc 2 [l5^ ,4 ft 2 7 + tofPHri]) 
+^c G l ([l8^ 2 7 ' + 18^' 2 H 2 A] + fc [6^' 2 HZ] + fc 2 [V«7 + ^ 2 r?]) , 
C2 i + S £ ^7 - 2^' 2 7' - 2^' 3 A) 



C4 1 



^ (fc [-12^ 3 H 7 ' - W 4 nA + 1 V 3 H 2 7 ] + fc 2 [<p' 4 a - V ' 4 Z] ) 
J_ (k [~lh^U 2 i - 15v' 5 H 2 A + 15^' 4 H 3 7 ] + 2fc 2 [-<p' 6 HZ + cp' 5 Ha]) 



pi 



c G — ( fc [-V«7' ~ 4v3 ,2 -HA + 6ip"H 2 -f] + ^fc 2 [<^' 2 ct - tp' 2 Z] 



M 3 

n G = -^i(fc [-^V + 3<p' 4 Ha - + fc 2 [4</ 3 H 7 - 6v?'V' 2 7 + ¥?' 4 ^ - <A 

+ ^ H/^'a + 12 V ' 5 H 2 a - l^'^Ua - 3</V' V] 

+ fc 2 [-12ip"(pf !t H'y + 12v3 ,4 -H 2 7 - 3^' 4 -H'7 + fy"^' 4 A + 6pVV - 6v?' 5 -H0] ) 
+^c G l (fc [-2^'VV] - 2fc 2 K 7 + ^]) . 



(45) 



(46) 



(47) 



I 

Note that the spatial derivative of the isotropic pressure p in not write it here. Finally, in fc-space, the perturbed Galileon 
fc-space is not needed in the CAMB code, which is why we do field equation of motion, Eq. d36b . reads 



= % (fc h" + 2iU + iff A! + ip'HA + 2<p"A] + fcV Z + fc 3 7 ) 

+ 7l3 ~? ( k [ 12 7'V'H + 12ip' 2 UA' - 18ip' 2 U 2 A + 36tp"tp"HA + 12tp"Hj' + l2y"U'i + 18^ 2 -H'A] 



M 3 a 5 



fc 2 [&ip' 2 UZ + 2^' 2 Z' + itp"(p'Z] + fc 3 [Atp'U'y - 2ip' 2 A + V'7]) 



+ ( k [ 5 V 2 %V - 108^ /2 -H 3 7 ' + 54(^' 3 H 2 ^' - I98ip' 3 1-L 3 A + 216ip"<p' 2 H 2 A + lOW'ftV 



c 5 1 



+108ip' 2 HH'~/' + l^^UU'A] + k 2 [-6(p' 3 H 2 Z + m<p\' 2 UZ + \2ip l3 H' Z + 12ip' 3 HZ'] 
k 3 [-lO<p' 2 H 2 j - l2tp' 3 HA - 4ip' 3 Hr) + 24ip"<p"H~f + l2ip' 2 H'~f + 



+^ -s (fc [-24Cy 3 H 4 7 ' - M5<p' 4 H 4 A + 60^' 3 H 3 7 " + 60ip' 4 H 3 A' + 300<p"<p' 3 U 3 A + l8G V \' 2 U 3 i 

+180<^' 3 'H 2 'H'7' + 225v?' 4 -H 2 -H'A] + fc 2 [-45^' 4 -H 3 Z + Wip"ip' 3 H 2 Z + 15(p' 4 H 2 Z' + mtp^UU' Z] 
fc 3 [-36<^' 3 -H 3 7 - l2tp' 4 H 2 7] - 15(p' 4 H 2 A + i^U'r, + 36<p"<p' 2 H 2 j + 24<p' 3 HH'~f + I2ip" ip' s Hr)] ) 

+ m^S ( k \- m2 ^" + 6 ^ n2A ' - Wh 3 a + i2tp"n 2 A + unnw + 2a v 'wh'a] 

fc 2 [6^'-H 2 Z + 4<p"HZ + Atp'HZ' + 4<p'H'Z] + fc 3 [2-H 2 7 - 4<p'HA + AH'-/ + 2y"rj\ ) . (48) 
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As another consistency test, we have checked that the con- 
servation Eqs. d2Tll22l i in fc-space, 

X' + (kZ - 3HA)(p + p) + m{ X + X p ) + kq = 0,(49) 
q' + AUq +(p+p)kA- k x p + \kU = 0,(50) 

are satisfied by the fc-space perturbed expressions derived 
above. 



TABLE I. The model parameters for the Galileon models studied in 
this paper. The C2 parameter is tuned to yield the required amount of 
dark energy today and its exact value depends on the choice of the 
initial Galileon energy density p v ^. 



Models 


C3 


C4 


CB 


CG 


Galileon 1 


12.8 


-1.7 


1.0 





Galileon 2 


6.239 


-2.159 


1.0 





Galileon 3 


5.73 


-1.2 


1.0 





Galileon 4 


5.73 


-1.2 


1.0 


-0.4 



1. Synchronous and Newtonian Gauge Equations 

Here, we present the recipe to write the CGI perturbation 
equations in the synchronous and in the newtonian gauge 13^1 . 

The perturbed Friedmann-Robertson-Walker line element 
in the synchronous gauge is written as 



1 fa\ 2 „ 



1 fa\ 2 „ 

2 (fc) ^ 



and the other CGI quantities are given by 



(55) 



ds 2 s 



a 2 (r) [dr 2 - (5ij + hf^dx'dx 3 ] . (51) 



Latin indices run over 1, 2 and 3, Sij is the delta function and 
the spatial perturbed metric /if- = /if- (x, r) is given by 



d 3 k e lkx 



kikjh s (k,T) 



-6 [ k^hj 



V s '(k,r) 



, (52) 



where a superscript 'S' denotes quantities in the synchronous 
gauge, x is the spatial position vector and fc = k/fc is the unit 
vector mode in the k-direction. The CGI and the synchronous 
gauge quantities are related by means of the following rela- 
tions 



A = 
r/ = 2$, 

H 

a = 0. 



($' - mi) , 



(56) 



We do not present the full perturbed field equations in the 
synchronous and Newtonian gauges because they are not used 
in our modified CAMB code. However, note that CAMB works 
in the cold-dark-matter frame where A = 0, which is equiv- 
alent to the synchronous gauge written in a slightly different 
formalism. 



IV. RESULTS 



A = 0, 
V = -277 s , 
_ h' s 
2k' 



(53) 



The line element in the Newtonian (also known as longitu- 
dinal) gauge is diagonal, described by two scalar potentials 
and $, and reads 



ds 2 N = a 2 {r) [(1 + 2*)dr 2 - (1 - 2<&)dx t dx l ] . (54) 

Written in this way, the perturbed line element is only appli- 
cable to the study of the scalar modes of the metric perturba- 
tions. The two potentials are related to the Weyl potential <f> 



In this section we present and discuss our results, which 
were obtained using a version of the CAMB code l34ll suitably 
modified by us to follow Galileon gravity models. 



A. Background 

We compute the evolution of the cosmological background 
using the Friedmann equation, Eq. d24b . and the background 
Galileon equation of motion given by taking the zeroth-order 
terms of Eq. d36t : 
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FIG. 1. (Color online) Evolution of the ratio of the Hubble expansion rates of the Galileon and ACDM models, H / Hacdm (H = 9/3), and 
of the Galileon field equation of state parameter w. The evolutions are shown for the four models of Table Ufor different initial conditions. 
In the Galileon 1, Galileon 2 and Galileon 3 panels, on the left-hand side from top to bottom and on the right-hand side from right to left, the 
lines correspond, respectively, to p v> i/p m ,i = {1CP 4 , 1CP 5 , 1CP 6 , 1CP 7 , 10 -8 }. The same for the Galileon 4 panels but for p v> i/pm,i = 

{i(r 4 ,i(r 5 ,5 x i(r 6 }. 



The value of the Galilean background energy density p v ^ at 
the starting redshift, which we take to be z, = 10 6 , is deter- 
mined through the zeroth-order part of Eq. ( |32"1 >. 

+ il^ 3 ] + I&cg & 2 e 2 ] , (58) 

by the initial values of the field time derivative ipi and the ex- 
pansion rate the latter being given by the fixed matter and 
radiation components via Eq. ( |23| ) (the Galileon background 
energy density should be negligible at early times). We spec- 
ify Oi using ttmo = 0.265 and f2 r o rj 8 x 10~ 5 for the present 
day values of the fractional energy density of matter and radi- 
ation, respectively Since we are assuming a spatially 
flat Universe we need the evolution of the Galileon field to be 
such that fi^o w 1 — fl m o ~ 0.735. This can be done by 



choosing appropriately the value of the C2 parameter by a trial 
and error approach. As a consistency test, we have checked 
that Eqs. d23l |25l ) are satisfied by the numerical solution we 
obtain from CAMB. Moreover, we have also checked that the 
background expansion solution from CAMB agrees very well 
with those in the literature l28l - [33ll and from an independent 
code written in Python by us. 

In this paper we focus on four different sets of Galileon pa- 
rameters which we list in Table U In l32tl (to which we refer 
the reader for further details on the background evolution of 
these models) it was shown that these choices of parameters 
are free of ghost and Laplace instabilities for initial condi- 
tions with Pip,i/ ' p m ,i ~ 10 -5 . Here we shall use this and 
other choices of initial conditions which have not shown any 
theoretical instabilities yielding therefore viable cosmological 
evolutions. In Table ITT1 we list all the initial conditions used in 
this paper with the values of the C2 parameter and age of the 
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TABLE II. The values of the parameter C2 and of the age of the Uni- 
verse for all the initial conditions used in this paper. The age for 
ACDM is 13.738 Gyr. 



P<p,i/ pm,i 


C2 


Age (Gyr) 


Galileon 1 


10 -4 


-27.00 


13.978 


1CT 5 


-27.49 


14.317 


10 -6 


-27.56 


14.366 


1(T 7 


-27.58 


14.374 


10 -8 


-27.51 


14.375 


Galileon 2 


10" 4 


-12.600 


13.614 


icn 5 


-12.846 


14.256 


5 x 1(T 6 


-12.857 


14.286 


KT 6 


-12.885 


14.357 


lO" 7 


-12.891 


14.372 


10 -8 


-12.892 


14.375 


Galileon 3 


10" 4 


-14.760 


13.854 




-15.122 


14.296 


1CT 6 


-15.179 


14.363 


lO" 7 


-15.188 


14.373 


10 -8 


-15.189 


14.375 


Galileon 4 


10" 4 


-14.186 


13.833 


10 -5 


-14.519 


14.285 


5 x HT 6 


-14.539 


14.312 




500 1000 1500 2000 

I 



FIG. 2. (Color online) CMB temperature power spectra for the 
Galileon 3 model with two different initial conditions and for ACDM 
(dashed black), together with the WMAP 7-year (squares) 0] and 
ACT (circles) l43"ll data. From top to bottom, at I = 500, the Galileon 
lines (solid) correspond to p v ,i/ ' p m ,i = {10~ 4 , 10~ 5 }, respectively. 

initial conditions does not propagate into the expansion rate 
which is only sensitive to changes in w for times sufficiently 
close to today when dark energy is non-negligible. 



Universe. 

Figure[T|shows the time evolution of the ratio of the Hubble 
expansion rates, H = 9/3, of the Galileon and ACDM models 
and of the Galileon field equation-of-state parameter, w = 
Pip I Pip, where 



Pip 



C4 

M<'< 



' A/a 
AJpi 



c 2 [m - 



CG 



L 



We* 



(59) 



is the background pressure (the zeroth-order part of Eq. d33l l. 

Figure [1] shows that, depending on the initial condition, the 
expansion rate can be faster or slower than in ACDM for dif- 
ferent times during the evolution. Another noteworthy as- 
pect of the background evolution is the possibility of having 
phantom-like dynamics, w < — 1, while avoiding the pres- 
ence of ghost degrees of freedom. This happens because of 
the Galilean shift symmetry which protects against these kind 
of instabilities. The initial values of p Vy i can have a great im- 
pact on the evolution of w: the lower p v ^ the more negative 
the values of w will be. The reason is that lower values of p v 
in the past will force the energy density to grow more dras- 
tically (w < —1) closer to today when the field starts to be 
driven towards the de Sitter attractor evolution l30l [32I r42tl 
(see I3ll , l33tl for expansion history observational constraints). 
However, for p v < 10 -5 , the strong dependence of w on the 



B. Linear perturbation results 

We now look at the physical predictions of the full lin- 
ear perturbation equations derived in the previous sections. 
We always use the best fit parameters from the WMAP 7- 
year data results n m0 = 0.265, n s = 0.963, H = 
100/ikm/s/Mpc (h = 0.71), fl k = 0, where n s and fl k 
are the spectral index and the fractional energy density asso- 
ciated with the spatial curvature. These values are obtained 
for a ACDM model but may be modified once a Galileon 
gravity cosmology is assumed. However, for the purposes 
of our analysis of linear perturbations, it is sufficient to con- 
sider these values and we will provide a revised fit of the 
WMAP 7 -year data in Galileon cosmology in a future work. 
The amplitude of the primordial curvature perturbations is 
A^(fc ) = 2.43 x 10~ 9 at a pivot scale k = 0.002MPC' 1 . 
We set the initial conditions of the Galileon perturbation 7 
and its time derivative to be zero, and have checked that the 
evolution of 7 is insensitive to the exact initial values. 

As a consistency test of the results that follow, we checked 
that the perturbed quantities we obtain from CAMB satisfy the 
fc-space conservation equations, Eqs. 



1. CMB 

In Figure [2] we plot the CMB power spectrum for the 
Galileon 3 model and ACDM together with the WMAP 7- 
yearJH (squares) and ACT ll43ll (circles) data. Figure[3]is the 



10 




FIG. 3. (Color online) CMB power spectra for the four Galileon models for different initial conditions and ACDM, together with the WMAP 
7-year data (squares) 101 and ACT (circles) l43ll data. In the Galileon 1 and Galileon 3 panels, from top to bottom, at I = 10, the lines 
correspond, respectively, to p v .i/ Pm,i = {10 _4 (not visible), 10 -5 , 10 -6 , 10 -7 }, ACDM. The same for the Galileon 2 and Galileon 4 panels 

but for p v ,i/p m ,i = {l(T 4 (not visible), 1(T 5 , 5 x lCT e , *1(T 6 , 1CT 6 , 1(T 7 }, ACDM and p v ,i/p m ,i = {l(T 4 (not visible), 1(T 5 , 5 x 
1CT 6 }, ACDM, respectively. 



same as Figure[2]but for the four models of Table Q] with a log- 
scaled x-axis which highlights the low-l region. The effect 
of the Galileon field in the CMB power spectrum is mainly 
two-fold. 

Firstly, the modifications of the expansion rate can shift the 
positions of the CMB acoustic peaks. The value of the ini- 
tial condition has an impact on the background expansion rate 
and hence on the distance to the surface of last scattering, 
which translates into different positions for the peaks. For 
sufficiently small values of p v .i/ Pm,i S 10~ 5 (not plotted in 
Figure|2]since they are indistinguishable from the p v ,i/p m .i = 
10~ 5 case) the Galileon 3 curves have essentially the same 
peaks as a result of the almost identical expansion rate (c.f. 
Figure [TJ- The same applies for the other models Galileon 1, 
Galileon 2 and Galileon 4. 

Secondly, the late time evolution of the gravitational poten- 
tial can be also different from ACDM, resulting in a modified 



signal of the ISW effect on the largest angular scales (low I in 
Figure[3]l. For instance, the choice p v .i/ Pm,i = 10~ 4 is com- 
pletely ruled out for all the models shown, since the spectrum 
at low / is larger than the observational data by several orders 
of magnitude. In this case, the ISW effect is so pronounced 
that it dominates over the first acoustic peak and can also have 
an impact on the second one. 

Lowering the initial amount of dark energy helps to recon- 
cile the models with the data. However, for Galileon 4 there 
is still too much power on large scales. Note that this model 
differs from Galileon 3 by having a non-vanishing value of 
cq and it is impossible to keep lowering the initial Galileon 
density {p v .i/Pm,i ~ 5 x 10 -6 ) as ghost instabilities start to 
appear. This may be a hint that the strength of the derivative 
coupling cg can have a crucial impact on the predictions. For 
all the other models (Galileon 1 to Galileon 3), for sufficiently 
small values of p V: i/ p m ,i, the dependence on the initial con- 
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FIG. 4. (Color online) Time evolution of the Weyl gravitational potential <j> for the four Galileon models and ACDM (dashed) for 
k = {1.0, 0.1, 0.01 and 0.001} /iMpc -1 . All the models have the initial condition p 9t i/p m ,i = 10" 5 . At a = 0.1, for the 
k — {1.0,0.1,0.01} /iMpc -1 panels, and at a = 0.4 for the k = 0.001 MVlpc -1 panel, from top to bottom the lines correspond, re- 
spectively, to ACDM, Galileon 1, Galileon 3, Galileon 2 and Galileon 4. In the k — 1.0 /iMpc -1 panel, Galileon 2 is the model with the 
oscillating evolution of <f) near today. 



ditions become less pronounced and the fit to the CMB im- 
proves. There are still differences from the best fit ACDM 
model and from the data at low I but since the errorbars are 
also larger due to cosmic variance, Galileon 1 and Galileon 3, 
e.g., are still compatible with the observations. 

It is interesting to note that the CMB power spectrum for 
the Galileon 1 and Galileon 3 models can be quite similar al- 
though their C3 and c.4 parameters are different. This shows 
that there are, to some extent, degeneracies in the Galileon 
model parameter space. On the other hand, changing only 
one of the Galileon parameters can also change considerably 
the CMB predictions. For instance, in the top-right panel we 
plot the CMB power spectrum of a model sharing all the pa- 
rameters of Galileon 2 in Table Q] except that C4 = —1.659, 
for pipi/pmi = 10~ 6 (dashed red). Note that c 2 also differs 
because it is tuned to yield the required amount of dark energy 
today, giving c 2 = 14.362. We see that by changing only C4 



the predicted CMB spectrum gets considerably closer to the 
data for the lowest values of I. It is also interesting to note 
that all the models have the value of C5 fixed and we expect 
a richer phenomenology if we allow this parameter to vary as 
well. 

To further understand the CMB predictions of the Galileon 
model at low I, we plot in Figure [4] the time evolution of 
the Weyl potential, <j>, which is the relevant quantity for the 
ISW effect. We show the evolution for different values of k. 
The variety of evolution can be very rich within the parame- 
ter space of the Galileon model and depends on the scale un- 
der consideration. The evolution of <f> agrees, to some extent, 
with the ACDM model during the radiation dominated era. 
However, in the matter era, while <j> is constant in the ACDM 
model, that is not the case for Galileon gravity and the grav- 
itational potential does evolve with time. In particular, we 
note a very pronounced variation with time of <fi for Galileon 
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4 during the matter era and today which explains why there is 
so much power at low I in this model (c.f. Figure [3}. More- 
over, for the models shown, the gravitational potential suffers 
an overall deepening with time 114411 . in clear contrast with the 
ACDM model where the gravitational potential gets shallower 
with the onset of the accelerated expansion. For some scales, 
the gravitational potential oscillates with time which might be 
reason behind the oscillatory signature in the ISW region of 
the Galileon models CMB power spectrum in Figure [3] 



2. Weak lensing power spectrum 

The weak lensing signal of the CMB anisotropies is deter- 
mined by the lensing potential tp, which is an effective poten- 
tial obtained by integrating the Weyl potential, <j>, from today 
to the time of last scattering ll45ll (see also ll46ll for a concise 
description and application to modified gravity theories). 

The angular power spectrum of ip is plotted in Figure [5] for 
the four Galileon models and we see that it can be noticeably 
larger than the ACDM result on all scales, as a consequence of 
the pronounced time variation of <f> in these models (c.f. Fig- 
ure|4]i. The Galileon 4 model is the one where the gravitational 
potential deepens the most with time and it is therefore the 
model with the most lensing power. The initial conditions also 
have an impact on the result, especially for p Vt i/ p m> i > 10~ 6 . 
For instance, for the case p Vl i/p m ,i = 10~ 4 (which is not 
plotted) the power is higher by several orders of magnitude 
for all the models. 

This is an important result and it shows that weak lens- 
ing measurements have the capability to place strong con- 
straints on the Galileon gravity model. In particular, the 
Galileon 1 and Galileon 3 models, which have CMB tempera- 
ture power spectrum predictions similar to that of ACDM for 
Pip,i/Pm.i = 10~ 6 (red line), nevertheless have very distinc- 
tive predictions for the power spectrum of the lensing poten- 
tial. 



3. Matter power spectrum 

Figure [6] shows the linear matter power spectrum predicted 
in the different models. We have chosen to plot the power 
spectra at redshift = 0.31, which is the median redshift of 
luminous red galaxies (LRGs) in DR7 from the Sloan Digi- 
tal Sky Survey (SDDS) H3|. A recent estimate of the power 
spectrum of LRGs is shown by the points with errorbars re- 
produced in each panel |T[]. By plotting the matter power 
spectrum at the same redshift as the measurement, there is 
no need to make any adjustment for the growth factor to com- 
pare theory to observation. However, since we are plotting the 
prediction of linear perturbation theory in real space, there are 
three effects which could be responsible for any discrepan- 
cies between the theoretical spectra and the measurement: 1) 
Galaxy bias. This is generally modelled as a constant shift in 
the amplitude of the power spectrum on large scales, though 
simulations show that the bias is scale dependent, particularly 



for highly clustered objects JUt]. 2) Redshift-space distor- 
tions. Using peculiar velocities to infer the radial distance to 
a galaxy introduces a systematic shift in the clustering am- 
plitude. Again, this can be scale dependent l49ll . 3) Non- 
linear effects. This includes the familiar mode coupling be- 
tween fluctuations on different scales, but also, in the case of 
the Galileon models, possible screening effects which could 
introduce scale dependent departures from the linear pertur- 
bation theory predictions. 

There are different lines of evidence which point to LRGs 
being biased tracers of the dark matter distribution. Interpreta- 
tions of the measured clustering of LRGs in terms of empirical 
halo occupation distribution models suggest that these galax- 
ies reside in massive dark matter haloes, with an effective host 
halo mass of « 1O 14 M HoHH]. At the median redshift of 
the LRGs, this suggests a linear bias factor of b ~ 2. Mea- 
surements of the three point correlation function of LRGs can 
be used to infer their bias, and also return b ~ 2 ll53l Hp] , 
For such a high bias, the amplitude boost from redshift dis- 
tortions on large scales is expected to be modest. LRGs are 
therefore expected to have a clustering amplitude that is ap- 
proximately four times higher than that of the dark matter on 
large scales. The measured power spectrum plotted in Figure 
[6] is an estimate of the power spectrum of the haloes which 
host LRGs, and is not directly comparable with the estimates 
of the LRG bias factor outlined above. Reid et al. |Q]] pro- 
cessed the LRG density field by "collapsing" LRGs in com- 
mon dark matter haloes, to reduce the small-scale "fingers of 
God" redshift space distortion. Hence, massive haloes which 
host more than one LRG are given the same weight as a halo 
which hosts one LRG. Therefore, the effective bias of a sam- 
ple of haloes weighted in this way will be smaller than the 
effective bias when retaining the weighting of the number of 
LRGs observed. If we compare the ACDM power spectrum 
to the measurement in Figure [6] we see that the effective bias 
of this sample is closer to b ~ 

Nevertheless, despite this complication, it seems reasonable 
to demand that in a viable model, the observed power spec- 
trum of LRG host haloes should have a higher amplitude than 
the linear theory matter power spectrum. This simple require- 
ment puts many of the Galileon model power spectra plotted 
in Figure [6] at odds with the observed power spectrum. For 
these models, the success of the comparison with the data de- 
pends sensitively on the value of p Vl i/p m ,i- For instance, the 
initial condition p v .i/Pm,i = 10~ 4 has an excess of power 
clearly incompatible with the observations, as it would require 
a bias parameter fj < 1. Lowering p v ^/ ' p m ,i allows a better 
match to the observations to be obtained and the results be- 
come less sensitive to the initial conditions (lower initial con- 
ditions have nearly the same prediction as p v .i/ p m ,i — 10~ 6 ). 
However, all the models still produce an excess of power when 
compared to ACDM indicating that the formation of linear 
structure can be highly enhanced by the modifications of grav- 
ity in the Galileon model, a conclusion in agreement with pre- 
vious linear perturbation studies in the literature l32i l42i PPill . 
The Galileon 4 model is the one with the worst fit, even for 
the lowest initial condition p v ^j ' p m .i — 5 X 10 -6 (recall that 
in this model lower initial conditions lead to the appearance of 
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FIG. 5. (Color online) Angular power spectrum of the weak lensing potential ip for the four Galileon models with different initial conditions and 
ACDM (dashed). In the Galileon 1, Galileon 2 and Galileon 3 panels, from top to bottom, the lines correspond, respectively, to p v ,i/p m ,i = 
{10 -5 , 10 -6 , 10~ 7 } and ACDM (the two smallest initial conditions are nearly indistinguishable in the Galileon 1 and Galileon 3 panels). The 
same for the Galileon 4 panel but for p v ,i/ p m ,i = {10~ 5 , 5 x 10~ 6 } and ACDM. 



ghosts). This indicates once again that the cq parameter can 
have a critical impact on the results. For k < 0.05 ZiMpc -1 , 
all the other models would agree very well with the data if 
6=1 and p v ^/ ' p m ,i 10 -5 . However, considering b > 1 
will increase the power on all scales, which could in principle 
be used to place strong observational contraints on Galileon 
models. 

We should stress, however, that when comparing the 
Galileon model with ACDM and clustering data one should 
be cautious about the validity of linear perturbation theory 
since the scale at which the Vainshtein screening effect be- 
comes important is not well known. For example, numerical 
simulations have shown that in other modified gravity mod- 
els such as the f(R) and dilaton f55l - l57ll . linear perturbation 
theory can fail even on scales as large as k ~ 0.01 /iMpc -1 
because of the screening ITlll57ll . As a result, a detailed study 
of the effects of the Vainshtein screening in Galileon gravity is 
necessary for a more complete comparison of the theory pre- 



dictions against the observations. This is beyond the scope of 
the present paper and will be left for future work. 



4. Clustering of the Galileon field 

We now turn the attention to the time evolution of the den- 
sity contrast S = p/p — 1. This is plotted in Figure [7] where 
we see that the Galileon density contrast (dashed line) can be 
large, being comparable with the dark matter (solid) and bary- 
onic matter (dotted) density contrasts throughout most of the 
evolution. This happens for all the scales considered including 
small scales such as k = 1.0 /iMpc -1 . 

This strong clustering of the Galileon field has a large im- 
pact on the evolution of the Weyl gravitational potential <fi 
which directly determines many observables such as the ISW 
effect (c.f. Figure[3]l, weak lensing (c.f. Figure|5]l and cluster- 
ing of matter (c.f. Figure|6]l. 
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FIG. 6. (Color online) Matter power spectrum at redshift zlrg = 0.31 for the four Galileon models with different initial conditions and 
ACDM (dashed), together with the SDSS-DR7 LRG host halo power spectrum QJ]. zlrg is the mean redshift of the LRG sample. In 
the Galileon 1 and Galileon 3 panels, from top to bottom, the lines correspond, respectively, to p v .i/p m ,i = {10 -4 , 10 -5 , 10 -6 } and 
ACDM. The same for the Galileon 2 and Galileon 4 panels panels but for p v ,i/p m ,i = {1(T 4 , 1CT 5 ,5 x 1CT 6 , 1(T 6 }, ACDM and 
p v ,i/p m ,i = {l(T 4 ,l(n 5 ,5 x 1(T 6 }, ACDM, respectively. 



One can also note that the Galileon density contrast starts 
to decrease with time close to the present day. This could be 
due to the rapid growth of the Galileon background density 
at those times {w < —1) which leads to a decrease of S v = 
P v /p v - 1. 



5. Quasi-static limit approximation 

In Figure [8] we plot the time evolution of the /c-space 
Galileon perturbation, 7 (dashed), along with the correspond- 
ing solution obtained in the quasi-static limit (solid). The 
quasi-static limit is the limit in which the spatial derivatives 
of the field are dominant over the time derivative ones. Prac- 
tically, this means neglecting all terms in the field equations 
that are suppressed by Ti 2 /k 2 or ip' /k 2 . 

As for the evolution of the density contrast and cf>, here there 



is also a strong scale dependence. Moreover, we see that even 
for near-horizon scales such as k = 0.001 /iMpc -1 the quasi- 
static limit can be a good (though not perfect) approximation 
to the full solution. In particular, in the Galileon 2 curves with 
k = 0.01 /iMpc - , one can see that the quasi-static approxi- 
mation agrees quite well with the full solution despite the os- 
cillations in the latter. The same model for k = 1.0 /iMpc -1 
has very rapid and pronounced oscillations, and therefore, the 
quasi-static limit cannot track the full solution because the 
time derivatives are not negligible. However, in this particular 
case, because the oscillations are rapid enough, only the aver- 
age value of the oscillating evolution is important and hence 
the quasi-static limit is still a good approximation. 

The quasi-static limit appears therefore to be valid for many 
cases in the Galileon model, especially when one is interested 
in subhorizon scales. However, it breaks down on superhori- 
zon scales and can lead to inaccurate predictions of the ISW 
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FIG. 7. (Color online) Time evolution of the dark matter, baryonic matter and Galileon density contrasts, S — p/p — 1, shown by the solid, 
dotted and dashed lines, respectively, for the four Galileon models for k = {1.0, 0.1, 0.01 and 0.001} hMpc~ . All the models have initial 
condition p v ,i/ p m ,i = 10 . All the models behave more or less in the same way for all the scales. In the k = 1.0 /iMpc -1 panel, Galileon 
2 is the model with the oscillating evolution of S v near today. 



effect and weak lensing signals. Moreover, as we can see from 
the lower-right panel of Figure[8] on near-horizon scales with 
k = 0.001/iMpc -1 , the error of this approximation can be 
a few percent, which is much larger than the numerical error 
of the CAMB code (which is at sub-percent level). For these 
reasons, we prefer to use the full numerical solution in the 
modified CAMB code. 



V. CONCLUSION 

We have studied the cosmology of Galileon gravity models 
at the linear perturbation level. For this we derived the full 
CGI perturbation equations using two independent methods: 
the normal procedure of linearising the full field equations and 
an alternative derivation that is particularly suitable for mod- 
els like Galileon gravity, where the shape of the Lagrangian 
is fixed by certain symmetries (e.g., there are no free func- 



tions such as the potential in quintessence and f(R) gravity 
models) and the field equations only contain up to second- 
order derivatives. The second derivation is particularly ap- 
pealing because it is much simpler than the first one, which 
is very lengthy and complicated for the full Galileon model. 
We checked that the two methods give the same set of per- 
turbation equations, and then solved these equations using a 
modified version of the CAMB code, which we tested by per- 
forming several successful consistency tests. 

Our code also solves the background expansion history in 
Galileon models and our results agree with those in the lit- 
erature. We find that the expansion rate in Galileon cosmol- 
ogy depends sensitively on the initial value of the Galilean 
energy density, especially if the latter is not small, e.g., if 
P<p,i/Pm,i ^ 10 -5 . Throughout the evolution, the expan- 
sion rate can be faster or slower than in ACDM and the 
Galileon equation-of-state parameter can cross the phantom 
line (w < —1) in a way which is free of ghost-like instabili- 
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FIG. 8. (Color online) Time evolution of the fc-space Galileon field perturbation 7 (dashed) along with the corresponding quasi-static limit 
(solid), for the four Galileon models for k = {1.0, 0.1, 0.01 and 0.001} hMpc -1 . All the models have initial condition p v ,i/p m ,i = 10 -5 . 
For the k — 1.0 MVIpc -1 panel at a — 0.3 and for the k = 0.1 ft.Mpc -1 panel at a = 1.0, from top to bottom, the lines correspond to the 
models Galileon 4, Galileon 2, Galileon 3 and Galileon 1, respectively. For the k — 0.01 /iMpc -1 panel at a = 1.0, from top to bottom the 
lines correspond to Galileon 2, Galileon 3, Galileon 1 and Galileon 4, respectively. For the k = 0.001 /iMpc -1 panel at a = 1.0, from top to 
bottom, the lines correspond to Galileon 2, Galileon 4, Galileon 3 and Galileon 1, respectively. In the k = 1.0 hMpc~ panel, Galileon 2 is 
the model with the oscillating evolution near today. 



ties. 

The modified background expansion translates into a differ- 
ent age of the Universe and distance to last scattering, which 
leads to a visible shift in the positions of the acoustic peaks of 
the CMB temperature power spectrum. The strongest effect of 
the Galileon field on the CMB temperature power spectrum, 
however, appears to be on the largest angular scales (low val- 
ues of I), where the full power receives a significant contri- 
bution from the integrated Sachs-Wolfe effect, which is due 
to the late-time evolution of the gravitational potential cf>. In- 
deed, we found that in Galileon models the gravitational po- 
tential evolves even during the matter dominated era and can 
undergo an overall deepening at late times. This is very differ- 
ent from the standard ACDM prediction that the gravitational 
potential is constant during matter domination and becomes 



shallower when the expansion of the Universe starts to accel- 
erate. The origin of the abnormal evolution in the gravitational 
potential can be traced back to the nonzero pressure perturba- 
tion and anisotropic stress of the Galileon field, which cause 
it to cluster strongly (comparable to the clustering of dark and 
baryonic matter species) on all scales. 

The evolution of the gravitational potential influences a 
number of cosmological observables, both directly and indi- 
rectly. In addition to the ISW effect, it also has strong impact 
on the growth of matter density perturbations (and therefore 
the linear and nonlinear matter power spectra), weak gravita- 
tional lensing and their cross correlations. In particular, we 
have shown that the Galileon model can predict considerably 
more power than ACDM for the weak lensing power spec- 
trum at all scales, even if their predictions for the CMB power 



spectrum more or less agree. Moreover, we have seen that in 
the low-Z region, the CMB power spectrum can oscillate with 
I, which can be attributed to the scale dependence in the evo- 
lution of the gravitational potential. Galileon models might 
also have distinctive predictions for the cross correlation of 
the ISW effect with the galaxy distribution. These are impor- 
tant observational signatures in the linear perturbation regime 
that can in principle help to distinguish the Galileon models 
from the standard ACDM paradigm. 

On the other hand, the sensitive dependence of the Galileon 
behaviour on the model parameters makes the phenomenol- 
ogy of the Galileon cosmology especially rich. For example, 
by tuning the parameters in the Galileon Lagrangian, one can 
get a CMB power spectrum which is very close to the ACDM 
prediction and therefore hard to distinguish by looking at very 
large scales. 

On subhorizon scales, we have seen that the linear growth 
of matter density perturbations can be significantly enhanced 
with respect to the ACDM results, even for those model pa- 
rameters that lead to similar CMB power spectrum. However, 
in Galileon models, the Vainshtein screening mechanism is at 
play and its potential influence on the clustering of matter is 
still to be properly understood. As an analogy, in other mod- 
ified gravity models such as the f(R) gravity, the chameleon 
screening effect has been shown to make the linear perturba- 
tion theory a poor approximation even on scales as large as 
k = 0.01 /iMpc -1 . We therefore conclude that a better un- 
derstanding of the true impact of the Vainshtein screening is 
necessary, before attempting a more rigorous confrontation of 
the predicted matter power spectrum with measurements of 
galaxy clustering. Such as study will be left for future work. 

Finally, we have seen that the quasi-static approximation 
for the evolution of the Galileon field perturbation serves as a 
good approximation on subhorizon scales for the models we 
have shown in this paper. It works reasonably well on near- 
horizon scales such as k = 0.001/iMpc -1 , with an error of the 
order of a few percent. However, for accuracy considerations 
we solve the full evolution equation of the Galileon perturba- 
tion in our code, which does not take much longer anyway. 

In conclusion, we have shown that the detailed study of the 
full perturbation equations unveils a rich phenomenology in 
Galileon gravity models. The full cosmological parameter 
space increases considerably in Galileon gravity. Even with 
current data, the indications are that strong constraints can be 
placed on this parameter space. In a future project we will use 
our modification to the CAMB software to carry out a formal 
study of the Galileon parameter space. 
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Appendix A: The Covariant Field Equations in Galileon Gravity tions are given by: 

The modified Einstein field equations and the Galileon field G>„ = K [T^ v + + T°f, + + T^ v + T^°] , (Al) 

equation of motion can be obtained by varying the action with 

respect to g^ v and cp, respectively. The Einstein field equa- where 



T, 



(A2) 
(A3) 



f'4 



9/j.u 



(CV) 2 Vav?V V - —R (V a ipV a cp) + 4n^V a 93V^^VaV^^ - AV x V a <pV x VfHpV a yN p ip 
-V x pV X (pV a V p pV a V p ip - R a pV a tpVP(pV X (pV X (p] 



-2V Q V ^V a V^ pV v (p + 8V X ipV P V x tpV p V {fl tpV „ } p - 2U(pV \<pV x <pV M V ' v <p 
~4\7 a \7 f! ip\7 a ^\7 f3 ip\7 tI \7 I ,ip-^RV X ipV x pV f ,pV u ip+^R f , u (y a (pV a tpf 



-2V ll \7 a <pV„V a pX7 x >p\7 x <p + IV^V^V^V^VVV^ + 2W t , a „p\7 a <p\7l 3 <p\7 x pV x <p] , 



(A4) 



T;i = j^9„u [0» 3 V A¥ >VV + 3(n<p) 2 V a V p (pV a <pVp(p - 3[tyVA¥>VVV«V j 8pV a V'V 
-6Dip\7 a \7 p ip\7 a \7 x p\7 pS7 X ip + 2V x ipV x pV a V^^V^V 7 ^V 7 V a ^ 

-3V«V /3 v5V"vVv p V (J v?V'VV'V + 6V q V /j <^V /3 V 7 <^V 7 Va^V q v5V a ^ - RDip (V x tpV x ipj 2 
+ ^i? a/3 V Q VV (Va^vV) 2 + 3Va^VVV>VVV 7 V>W Q7(3o 
+ ^ [(□v 3 ) 3 V Al ^V^ - 3(\3p) 2 V x pV x tpV fl V^ - 6(D</>) WVaV^V^ + 6npV x <pV x tpV fl V a <pV u V a ip 



-6n^v Q v / 3v?v Q >\7 / Vv M v I ,^ - 3n^v Q! v^v Q v , Vv AI < / 3V I ,</j + 6npV a pV p pV f _ t v a pV u Vfj(p 

+ 12n^V Q v3V Q V' 3 v5V /3 V (AI <y9V, /) ^ + 3VA^V A (^V Q V ;3 ^V Q V ,3 v5V A1 V t/ <y9 - 6VA</>VVV a V / VV M V a </>V„V / 3p 
+&V a V p <pV a tpV ptpX/ ^X/ x pX/ u V x <p + 6V Q V A ( / 3V /3 VA<^V Q ^V' 3 ( / 3VpV 1 ,( ( 5 + 2V a S7 p tpV /3 V x ipV x X7 a tp\7 ^ipV u tp 

-12V Q ^V a V A ^V T ^VAV (A1 ^V iy) V CT (p + !□¥> (V x <pV x <p) 2 - 3 (V x <pV x ip) 2 R^V v) V c <p 
+R (V x ipV x <p) 2 V M V„y> + 3n<pVx<pV x <pV*<pV p <pWr aV (i - W\<pV x <pV*<pV p V y <pVfr<pW v ) PBn 



-6VA^V A ^V a ^V 3 ^V 7 V (A1 ^W, )Q(37 - 6Va^VVV q ^V q V^V 7 ^W ( f 7 



(A5) 



Mpi 



= -^c G ((D^) 2 - V Q V^V a VV) + 2V„V A ¥>V„VV - 2D^V M V^ 



+2R x{ll V v) (pV x ip - R al3 V a pVf3<P9^ + 2W a ^ u V a tpV p tp - ^i?V^V^ + ^V^Va^, 



(A6) 



The Galileon field equation of motion is given by: 



o = c 2 n<p + 2-^ [(Dip) 2 - v°vVv a v^ - i? Q0 v>vV] 



M 6 



2{\3<p) 3 - 6QpV Q V (3 <^V a V / V + 4V Q V /3¥ A7 /3 V'VV 7 V> - -R57 a <pV%\J a \J [3 Lp 
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--ROipV\ipV x tp + AR a pV a V f3 ipVxipV X ip + 4i? a/3 V t VV>V /3 VAV? - 2D< / ?i? a/3 V>V / V 
o 



~ \(aip) 4 - 6(Q£) 2 V Q .V /3 y>V Q! V / V + 3 (V Q V /3 ^V Q V /3 <^) 2 - 6V Q VVV/jV„v5V"V A 
M J L 

+8D^V ct V /3 ^V, 3 V A <y9VAV Q ^ + 6aip\7 x ip\7 X ipR ap \7 a \7 l3 ip - 2i?(Q y j) 2 V A ^VV 
1 



+ ^RR a[ 3 V a cpV P (pV x cpVxV> - 6i? Q(3 V A <y9 V A f V° V 7 99 V 3 V CT </? + 2R V A 99 V a (p V Q V 9 99 V Q V p cp 

~R p(7 R pa (Va^V a v?) 2 + \R 2 (Va<^V A (^) 2 - ii?i? Q/3 V Q (^V' 3 ^V A ^VA^ 

+6DipW pac7 pV p V a (pW a ipW ip + l2W pa p„\' p V' 7 ipV l3 V x ipV a ip\'x<P 
+3W pct(3ff V' , V>V Q V^VvVVA< / J + 6W Q(9(T/3 VVV' t ^V q V a v3V /3 Vav 

+2^|c G G Q/3 V Q VV (A7) 



The usual equations presented in the literature f lfl3l l32ll . 
e.g.) are related to ours via the following Riemann tensor ex- 
pansion 



Rp.vuf3 — ^ {9p.aRuf3 + 9upR p a ~ 9^Rua ~ 9vcxR p fi) 

+W pi y a p - -R {g^ a g v p - 9n[i9va) , (A8) 
which cancels some of the terms originally derived in loll . 



Appendix B: Alternative Derivation of the Perturbed Equations 

In this appendix we present an alternative derivation of the 
perturbed equations in Galileon gravity. This method requires 
only the knowledge of the Galileon equation of motion and the 
assumption that all the field equations do not contain deriva- 
tives higher than second order, the latter being satisfied by the 
theory of Galileon gravity by definition. 

If the above requisites are satisfied, then it is easier to derive 
the perturbed components of the Galileon energy-momentum 
tensor using the new method rather than from the complicated 
Galileon Lagrangian. In the latter case, one has to first de- 
rive the full expressions of the energy-momentum tensor (see 
Appendix [A), which itself could be a considerable amount of 
work. 

The spirit of this derivation follows the general method in- 
troduced in llslll and generalised later by l!59Tl6lll . However, 
here we work within the framework of covariant and gauge- 
invariant perturbations, and consequently the mathematical 
description looks different from those works. 

To lighten the notation, in this appendix we neglect the su- 
perscript G in the dynamical quantities for the Galileon field. 



1. The Method 

As we have seen above, the quantities p, p, and ir^ have 
contributions from both normal matter and the Galileon field 
ip. Here, let us first look at the most general forms that p, p, q p 
and ir^ for the Galileon field can take. The arguments are as 
follows: 

1 . Eq. ( fT8l ) contains time derivatives up to first order (in 8) 
and spatial derivatives up to second order (in R). If we 
want to keep this property, the Galileon energy density 
p can contain 8, R, tp and Dip, but not their time deriva- 
tives or gradients. It cannot contain quantities such as 

\7 v cr piy , which involve higher order derivatives. If it 
contains V ■ A, then according to Eq. ( 1211 p must con- 
tain (V - A)' which involves third-order derivative and 
hence it is not allowed. 

2. According to Eq. ( |2TI ). the Galileon pressure p can con- 
tain 8, 8, tp, R, (p. It can also contain V • A without 
changing the structure of Eq. dT3T >. Quantities such as 

R, (V • A)' and \7 P \7 V a ^ are not allowed as they con- 
tain higher-order derivatives. 

3. The Galileon field peculiar velocity q p can contain 
V M #, ^J v a pv , \/ v vo pv , VjuV 3 and \7 p ip, but not then- 
time and spatial derivatives. If it contains A„, then 
Eq. ( l22l cannot hold without involving derivatives 
higher than second order. It cannot contain V„ii be- 
cause this has third order derivatives. Cp and 8 are not 
allowed because otherwise either p or ir^ has to con- 
tain third order time derivatives, according to Eq. ( l22l . 

4. The Galileon anisotropic stress tensor n pu can contain 
cr^v, 6-^,8^, V^V^ and W^A^, but not their time 
and spatial derivatives. 

Based on the above analysis, we can decide which terms 
can appear in the expressions of p, p, q p , ir^ for the Galileon 
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field. More explicitly, up to first order in linear perturbations, 
we have 



B p ip a 6 b - 2 R + C p cp a - 1 b - 1 ncp, 



p = A p tp' u T u p 

■ p a-l e b-l + Bpt pa,Q 9 b-2 + c ^ 



(Bl) 



p = Aptptp' 

p a - 2 e h - 2 Uip + E p ip a ~ 1 90 b - 3 np 



9m 



+D pl p<p a 

+F p (p a ~ 1 b ~ 1 Dp + G p ipp a ~ 1 e b - 3 V ■ A 
+H p p a 99 b - i V ■ A + l p ip a 9 b - 2 V ■ A 
+J p <pcp a - 1 6 b - 3 R + K pl p a 06 b - i R 
+L p <p a 8 b - 2 R, 

A q ip a ~ 1 6 b ~ 1 V fl ip + B q tp a - 1 6 b V fll p 

D q p a e b - 2 v»<Tp„ 



+C q p a d b - 2 Vp6 + D ( 
+E ql p a 6 b - 2 V' J vjp V 
iTp„ = A^ a B b - 2 £p V + Bi T ip<p a - 1 6 b -%„ 
+C„<p a 99 b -%„ + D^ a d b - 2 bp V 

-wv a - 1 9 b - s & tiV + F^ a ee b ~ A bp V 



(B2) 



(B3) 



y CLQb— 1 



+GnP 



op 

|6-3 



(Tpv 



- J^p a B b ~ 2 V { pA v 



K 7r pp a - 1 8 b - 3 V { pA v} + L^ a ee b - i v {ll A v) 

M^ a - 2 6 b - 2 V { pV v) p + AT^-^V^V^ 

(B4) 



+o^ a ~ 1 ee b ~ 3 w { pV v) p, 



in which A P PtgtW , Bp.p.q.-n-, • • • are constant coefficients and 
a, b are dimensionless constant power indices. Note that to 
write down the above equations we have used the fact that 

1. all terms in the expressions must have the same mass 
I 



dimension and 

2. the power of tp (with <p, <p and Dip counted in) must be 
the same in all terms, 

which must be true if the dynamical quantities are to be de- 
rived from the Lagrangian densities C1—C5 that are specified 
in Eq. ©. 

When using the above expressions, we require that all terms 
must not contain negative powers of 9 (which will never ap- 
pear when varying £1-5 with respect to the metric gpv). For 
example, if b = 3, then K p should be set to zero. 



2. Application of the Method: the C4 Term 

Here we illustrate the application of our method for the par- 
ticular case of the a term. The Lagrangian £4 is sufficiently 
complicated to highlight how much simpler this method can 
be. For this term, we know from the background expression 
of the energy density (or equivalently the Galileon equation of 
motion) that a = 4, b = 2, and so we can write 



p = A p p 4 9 2 + B p ip 4 R + C> 3 60<p, 



(B5) 



A p p<p 3 6 + B p (p 4 6 + C p (p 4 6 2 + D p (p(p 2 Up 
+F p (p 3 6\3ip + I p p 4 V ■ A + L p p 4 R, (B6) 
A q (p 3 6S7p(p + B q (p 3 e 2 Vpp 
+C q p 4 Vp9 + D q p 4 V»(jp y + E q p 4 i] v vop V <&l) 
Kpu = A^Epy + D^bpv + G 7T (p 4 9ap V 
+H 7T (pp 3 (7p U + J T (yj 4 V \pA v ) 
+M 7T ^<p 2 V { pV u) tp + N^ 3 9V { pV v) ip. (B8) 

Substituting these into the conservation equations ( F2T1|22] |. we 
find 



{AA P + Ap)Cpp> A 9 A + (2A P + Bp)<p*60 + (A p + C p )tpW 3 + 



( \ b p + L v 



(3C P + Dp)Cpp 2 B + C p p 3 9 + ( -A g + B q + F p + G p ) <p 3a ' 



atp + (A q + c p )p 3 e[n<p 



R 



-Bp + C q 



p 2 ne- 



A q + I p ~ -Bp 



(A q - A p )<p 3 8 (V^) 



<p A 0V ■ A + (D q — 2B p )p 4 ij^\/ lJ (jp U + {E q - 2Bp)p 4 ij^ v Wp U = 0, (B9) 



>p 2 9 + {A q - 4B p )p 3 9 + ( B q + ^A q - 4C P - ^A p ) <p 3 ° 



3ptp 2 8 z + tp 3 3 + 2<p d 69 Vpip + {A n - \2Lp)^\I v Ep V + {D q + D n — 6L p )vj 4 Vct 



+B q 

-{E q - 6L p - J n )ip 4 (Vw F 



> 3 + ( \o q + Id„ + G^ - AL p ) ^ 



V V Gp V 



{AE q + AU)pp 3 + [ -E q + N„ - ALp - J« ) p 4 



Vmpv + {A p -B q + Cp)ip 4 2 Ap 



iC q -Bp)p 4 (VpO)' 



( ^C q - 2C P - ±B P ) p 4 < 



V p 
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I p — ALp — g^TT 



D p - -M n ) <ptp* 



V^Qp = O.(BIO) 



From the background expression of p and p (or equivalently 
the background Galileon equation of motion together with the 
energy conservation equation) for the ca term, we find 

A p = ^A, A p = -4A, Bp = -A, C p = (Bll) 

in which A = C4/M 6 . This can be done by equating the first 
three terms of Eq. dB9l ) to the background Galileon equation 
of motion 



A v = 2D V + 3A, 
1 



G, 



- A, 

1 



3 D * X 

2 2 
F n = --D w --X. 



(B13) 



3^<9 2 



2(p 2 6 



2/)3 



= 0, 



(B12) 



which can also be used to eliminate the terms containing W^ip 
in Eq. dBlOt . 

Because we have already used the Galileon equation of mo- 
tion in Eq. ( IB 10b , the remaining terms on the left-hand side of 
this equation must cancel amongst themselves. In addition, 
for Eq. (|B9t to not contain higher-order derivatives, we must 

set the coefficients of (Oip^j , 06, V^VV^ and W^^zu^ 

to zero. This gives us 



B 



and 



C p = -Dp = —A q = 4A, C q = -Ip = -A, 
, - 2A, H n = M n = -6A, D q = E q = 2B P = |a, 



L"n — ~D~ 



Unfortunately, some coefficients cannot be fixed unambigu- 
ously, and here we have expressed all those coefficients in 
terms of D^. This indicates that perhaps the Galileon model is 
not the only one which gives the perturbed energy-momentum 
tensor as in Eqs. dBUB2IB3IB4b . To solve this problem, we 
can use the perturbed Galileon equation of motion to fix the 
free parameter. Of course, this does not necessarily mean that 
we have to write down the full perturbed equation of motion. 
Indeed, we only need to know the coefficient L p or F p . 

The Galileon equation of motion can be read from the re- 
maining terms of Eq. (|B9l l, from which we find that the ratio 
of the coefficients of (p(p 3 d 2 and ip 4 6R is 1/12 + D X /(36A). 
On the other hand, the value of this ratio can also be easily 
calculated by explicitly perturbing the Galileon equation of 
motion where we find it to be 1/18. As a result, D^ = — A 
and all the coefficients are now fixed. The components of the 
energy-momentum tensor of the C4 term are: 



P 
P 



C4 

M 6 

CA 

M 6 

CA 

M 6 

CA 

M 6 



^e 2 +'^ 4 R + ^ 3 eOif 

-ACp^o - (p 4 e - \ip A e 2 - 4<^3 2 n< / 3 - -ifeoip + ^ 4 v • a + -^^r 

2 9 12 



-40 3 6>V^ + 2<p 3 e 2 V li <p - £ 4 V M + -<ys 4 V<>^ + Wp U ) 

4 2 



(B14) 
(B15) 
(B16) 
(B17) 



and the fully perturbed Galileon equation of motion becomes 

= 6fiip 2 8 2 + Aip 3 m + 2tp 3 9 3 - 4ip 3 6V ■ A + 



30 + fiQ 



<f 2 R- 



26 

8(fiLf>e + 4ip 2 e + —ip 2 e 2 

9 



n<p, 



(bis) 



which is in agreement with the C4 terms in Eqs. d3"2l-l35Tl and We have applied the same method to all other terms, and for 
Eq. (l36l l. respectively. all of them the resulted equations agree with Eqs. (f32]-[35]l and 
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Eq. (l36fl Note that in this new method the different terms of 
the Galileon field can be worked out in a unified way, which 



further reduces the computational effort. With certain modi- 
fications, the method should be applicable to the generalised 
Galileon model fill as well. 
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